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ABSTRACT 

Banerji et al. (2009) suggested that top-heavy stellar Initial Mass Functions 
(IMFs) in galaxies may arise when the interstellar physical conditions inhibit low- 
mass star formation, and they determined the physical conditions under which 
this suppression may or may not occur. In this work, we explore the sensitivity 
of the chemistry of interstellar gas under a wide range of conditions. We use 
these results to predict the relative velocity-integrated antenna temperatures of 
the CO rotational spectrum for several models of high redshift active galaxies 
which may produce both top-heavy and unbiased IMFs. We find that while ac- 
tive galaxies with solar metallicity (and top-heavy IMFs) produce higher antenna 
temperatures than those with sub-solar metallicity (and unbiased IMFs) the ac- 
tual rotational distribution is similar. The high-J to peak CO ratio however may 
be used to roughly infer the metallicity of a galaxy provided we know whether it 
is active or quiescent. The metallicity strongly influences the shape of the IMF. 
High order CO transitions are also found to provide a good diagnostic for high 
far-UV intensity and low metallicity counterparts of Milky Way type systems 
both of which show some evidence for having top-heavy IMFs. We also compute 
the relative abundances of molecules known to be effective tracers of high density 
gas in these galaxy models. We find that the molecules CO and CS may be used 
to distinguish between solar and sub-solar metallicity in active galaxies at high 
redshift whereas HCN, HNC and CN are found to be relatively insensitive to the 
IMF shape at the large visual magnitudes typically associated with extragalactic 
sources. 

Subject headings: astrochemistry-galaxies: active - ISM: abundances - ISM: molecules 
- stars:formation 
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Introduction 



In iBanerji et al.l (120091 ) (B09 hereafter) we investigated timescales relevant for low-mass 
star formation for galaxies with a variety of different physical conditions. That study allowed 
us to infer in a crude way the physical properties of galaxies in which low-mass star formation 
is likely to be impeded possibly resulting in a top-heavy Initial Mass Function (IMF). We 
concluded that low-mass star formation is likely to be inhibited in low-metallicity systems as 
well as systems with high cosmic ray ionisation rates and FUV radiation field strengths where 
the magnetic pressure dominates over the thermal pressure. Such systems may therefore have 
a high-mass biased IMF. We also concluded that active galaxies at high redshift with sub- 
solar metallicities show enhanced low-mass star formation resulting in an unbiased IMF. 
However, similarly active galaxies with solar metallicities may have low-mass star formation 
impeded by magnetic pressure resulting in a high- mass biased IMF. 

The aim of this follow-up paper is to now consider the interstellar chemistry associated 
with extragalactic regions with a variety of different physical conditions. By looking at the 
chemistry associated with the different extragalactic systems studied in B09, we can identify 
possible molecular tracers of different types of IMF that can be used to direct observations 
with future projects such as ALMA. In particular, we consider the rotational transitions of 
^^CO as well as molecular abundances associated with different physical conditions. In § 
12] we describe our model for interstellar chemistry. In § E] we consider changes in the CO 
antenna temperature for the various rotational transitions for various physical conditions in 
the collapsing cloud. In § H] we look at changes in the molecular abundances for different 
physical conditions and as a function of depth within the cloud. Finally, we summarise our 
results in § O 



The Model 



All results in this paper are obtained b y modelling the gas-p hase chemistry inside a 
molecular cloud using the UCL_PDR code flBell et al.l |2006| . 120071 ) for photo n-dominated 



regio ns which has already been benchmarked against several other PDR models (IRollig et al. 



20071 ). The model is both time and depth dependent but has been run here at steady state 
corresponding to a time of IGyr. The cloud is treated as a semi-infinite ID slab illuminated 
by an unidirectional flux of far ultraviolet photons. The code solves simultaneously the 
chemistry, thermal balance and radiative transfer within the cloud. The chemical network 
considered here has 83 spe cies that interact in over 1300 reactions whose rates are taken 
from the UMIST database Jle Teuff et al.ll200oh . 
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In the models considered in this work, the hydrogen atom number density, n, metalhcity, 
^, cosmic ray ionisation rate, ( and FUV radiation field strength, x are the input parameters 
that are varied. The turbulent velocity of the cloud is fixed at 1.5 km s~^ and the dust to gas 
mass ratio is assumed to scale linearly with metallicity. All the models used in this paper 
have already been presented in B09. 

Central to much of the work presented in this paper are the model predictions for pa- 
rameters associated with the rotational lines of CO and we therefore describe how these are 
obtained in more detail. The gas in PDRs mainly cools through emissions from collision- 
ally excited atoms and molecules as well as through interactions with cooler dust grains. 
The code calculates the emission from the [01], [Oil] and [CII] fine structure lines as well 
as the CO rotational li nes at each depth a nd time step. This is done using the escape 
probability formalism of Ide Jong et al. to solve the radiative transfer equations. Non 



local thermodynamic equilibrium (non-LTE) level populations are determined by solving the 
equatio ns of statistical equihbrium using rate coefficients fo r CO collisions with H2, H and 



He from iFlower fc Launayl (Il985l ) , ICreen fc Thaddeud (Il976l ) and ICreen &: ChapmanI (Il978l ) 



respectively. The solutions of the radiative transfer equations give the emissivity. A, the 
intensity, I and the opacity, r for every line at every depth and time step. Several of these 
outputs are used throughout this paper. For further details of the UCL_PDR code and the 
models considered in this work, please refer to B09 and references therein. 

Although this paper makes references to models of extragalactic systems, we emphasise 
that we are modelling the various PDR components that make up external galaxies rather 
than the galaxy as a whole. Our chosen input parameters therefore mimic the average con- 
ditions in individual PDRs within a galaxy. We choose to consider outputs at two visual 
extinctions of ^ 3 and ~ 8. These roughly correspond to the translucent and dense 
PDR gas components within a galaxy. Atomic and molecular emissions from different com- 
ponents within the galaxy will contribute independently to the line emissions detected from 
the galaxy as a whole and studied in detail in this paper. Also, due to the fact that the 
PDR is modelled as a semi-infinite slab, a given visual extinction will correspond to different 
physical distances within the cloud under different physical conditions as the y4„ depends on 
both the number density and the metallicity of the cloud in our models. We come back to 
this point later in the paper. 



3. ^^CO Rotational Transitions 



Carbon monoxide is one of the most commonly used tracers of molecular gas and its var- 
ious rotational J transitions have been detected in many extragalactic systems both nearby 
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2003 



(iBradford et al.l 12003: iBayet et al 



Nesvadba et al 



2009 



200' 



Daddi et al 



2009 



200t) and in the distant Universe (IBertoldi et al. 



Knudsen et al.ll2009l ). In Milky Way type con- 



ditions, the column density of molecular hydrogen along a line of sight is often related to the 
velocity integrated antenna te mperature, for the CO(l-O) transition, using the CO-to-H2 
conversion factor or X factor (jStrong fc Mattoxlll996l : iDame et al.ll200ll ). The CO(l-O) line 
is used as it is known to be the dominant CO cooling line with the highest antenna tem- 
perature under physical conditions appropriate for a Milky Way cloud. However, this is not 
necessarily the case for high redshift systems with significantly d ifferent physical conditions 



and higher order transitions of CO may become more relevant (jSakamotd Il999l : iBell et al. 



20061 . 120071 ). It is alre ady known that galax ies with elevated cosmic ray fluxes have larger 
high- J CO line ratios (IMeijerink et al.ll2006l ) and that the higher order CO transitions trace 
denser gas than their lower order counterparts. In this section, we examine CO rotational 
line transitions and ratios in the context of the stellar initial mass function. We consider 
which of the CO transitions, if any, can be used as diagnostics for regions of both enhanced 
and inhibited low-mass star formation as studied in B09. We conduct this study both for 
systems where we vary one physical parameter at a time while keeping the others fixed at 
their local values as well as the three high redshift systems modelled in B09 whose input pa- 
rameters are summarised in Table [1] along with the IMFs associated with them and inferred 
in B09. 

The UCL_PDR code outputs the integrated line intensity of the CO rotational lines with 



1 < Jupper < 9 in units of erg s" 



-1 —2 

cm 



sr ^ after solving the radiative transfer equations. 



The intensity is related to the emissivity A by 



1 = ^1 mdL 



where L is the depth into the cloud and the factor of l/27r takes into account the fact that 
the photons only emerge from the outer edge of the cloud due to the cloud being modelled as 
a semi-infinite slab. Integrating over larger depths will therefore obviously result in higher 
intensities. 

In order to calculate the relative velocity integrated antenna temparature from the 
intensity the following equation is used: 



10' 



c^I 



Kkms ^pc ^ 



(2) 



where c is the speed of light and /c^ is the Boltzmann constant. As the UCL_PDR code 
models the PDR as a semi-infinite slab and the visual extinction, in the slab is related 
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to both the number density as well as the metallicity, a fixed Ay will correspond to different 
physical distances, L, for different values of the density and metallicity. In order to derive 
values for the relative temperatures that are comparable and independent of the size of the 
emitting region, we also divide the relative velocity integrated antenna temperatures by the 
distance in parsec corresponding to the at which the temperature is calculated, La^ ■ This 
distance will obviously be much larger for low density and low metallicity systems for a fixed 
Ay. Our figures therefore represent the relative velocity integrated antenna temperatures for 
the different ^^CO transitions per unit parsec except in Figure E] where we plot the relative 
velocity integrated temperatures function of A,,. 

Note that the numbers calculated using Eq [2] are important only in terms of observing 
trends. In order to match the absolute values of the theoretical velocity integrated antenna 
temperature to those from observations, one has to consider several other factors. Firstly, all 
models in B09 assume a turbulent velocity of 1.5 km s"^ typical for a giant molecular cloud in 
the Galaxy. In reality, observed lines of CO in extragalactic sources typically have widths of 
several 100 km s~^ due to contributions from several PDR regions within the galaxy and the 
temperature has to be scaled accordingly. Secondly, we need to account for a surface filling 
factor which takes into account the size of the source as well as the telescope beam. Finally, 
all observational results will be affected by factors such as the atmospheric conditions and 
the telescope efficiency. One way of comparing our theoretical predictions to observations 
however, is to consider line ratios rather than the intensities and brightness temperatures of 
individual lines. Assuming the emission from both lines comes from the same clouds, the 
various factors discussed above should cancel out when computing the ratio of either the 
intensity or the integrated temperature. We therefore compare our theoretical predictions 
to observed line ratios in § 13. 6[ 

In the following sections, we list in all cases our relative integrated temperatures per unit 
distance. We emphasise that the results presented here are only useful in terms of observing 
trends and should not be compared to absolute values obtained through observation as they 
do not take into account the factors already stated above. 



3.1. Density 



We vary the hydrogen atom number density of the molecular cloud between 100 cm~^ 
and 10^ cm~^ while keeping the metallicit y, cosmic ray ionisation rate and FUV fiux fixed 
at ,^ = ,^0, C=10~^^ s~^ and x=l-7 Draine (jDrainelll978l ). The effect on the relative velocity 
integrated antenna temperature, J {TAY'^''dv, of the different CO rotational transitions, is 
illustrated in Figured! The CO(l-O) line is seen to be the brightest line at all densities and 
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the CO SED falls off with increasing J. Higher density clouds also produce brighter antenna 
temperatures as expected. As the is increased, the relative integrated temperatures per 
unit parsec drop, particularly for the higher density clouds. This is telling us that in these 
high density clouds, most of the emissivity in these CO rotational lines is coming from the 
lower Ay regions. This is also the case for the high-J transitions in the lower density clouds. 

3.2. Metallicity 

We vary the metallicity between ^ = 1O~^^0 and ^ = 3^© while keeping the density, 
cosmic ray ionisation rate and FUV flux fixed at n=10^ cm~^, (^=10~^^ s~^ and x=l-7 
Draine. The dependence of / {TAY^^dv on the metallicity is illustrated in Figure [2l In BOO 
it was found that low-mass star formation is likely to be inhibited at very low metallicities 
resulting in a top-heavy IMF, provided the cloud is at moderately high densities and low 
cosmic ray fluxes and FUV field strengths. The form of the CO SED is quite different at 
different metallicities. While the higher metallicity clouds produce SEDs that peak at the 
lower rotational numbers and drop off sharply after that, the lower metallicity clouds have 
flatter SEDs. For the highest rotational transitions with Jupper > 7, the sub-solar metallicity 
clouds actually produce emission that is brighter than that from solar metallicity clouds. 
As we approach sub-solar metallicities, the ratio of the high-order CO transitions to the 
peak CO transition tends towards unity whereas at solar metallicities it is much smaller. 
The physical conditions in sub-solar metallicity counterparts of Milky Way type galaxies are 
thought to be appropriate for the production of top-heavy IMFs. 

3.3. Cosmic-Ray Ionisation Rate 

The cosmic ray ionisation rate is varied between C=10~^^ s~^ and (=10~^^ with the 
rest of the input parameters fixed at n=10^ cm~^, ^ = ^0 and x=l-7 Draine. The results 
are shown in Figure [31 

Apart from the highest cosmic-ray ionisation rates, the CO (1-0) line dominates the 
emission in almost all cases. The relative antenna temperature increases with increasing 
cosmic ray ionisation rate. At C = 10~^^ s~^, the CO (2-1) transition just starts to dominate 
over the CO(l-O) line but by C = 10~^^ s~^, the CO(5-4) line is now the brightest CO 
line. This line can therefore be used to trace high density systems at solar metallicity 
with extremely high cosmic-ray fluxes such as sources with AGN or a starburst nucleus. In 
B09 it was noted that systems such as this show some evidence for having top-heavy IMFs 
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depending on the exact core temperature and magnetic field within these galaxies. In some 
active galaxies with physical conditions similar to that considered here, the magnetic pressure 
created through ambipolar diffusion may be sufficient to halt core collapse and formation of 
low-mass stars. 



3.4. FUV Radiation Field Strength 

The FUV radiation field strength is varied between x=l-7 Draine and x=l-7xl0'^ Draine 
with the rest of the input parameters fixed at n=10^ cm~'^, ^ = ^0 and C=10~^^ s~^. The 
dependence of the CO SED on the FUV radiation field strength, is illustrated in Figure HI 

Increasing the FUV radiation field strength results in similar behaviour as increasing 
the cosmic ray ionisation rate as in B09. Higher-order lines become more and more dominant 
and the antenna temperature increases with x- At the maximum FUV field of ^ = 1.7 x 10^ 
considered here, it is the CO (8-7) line that is dominant but for a FUV radiation field that 
is an order of magnitude weaker, the CO (3-2) line dominates. The high-order CO transi- 
tions can therefore be used as tracers of extremely high FUV radiation fields associated with 
the presence of massive stars. Such systems may also have top-heavy IMFs depending on 
whether the magnetic pressure is sufficient to halt the collapse and formation of low-mass 
stars. 

Note that due to the large parameter space that we have tried to model, there will 
inevitably be degeneracies and the CO SEDs produced by models with x = 1700 and 
C, = 10~^^s~^ may be observationally indistinguishable. However, we have associated large 
values for both these parameters with active galaxies producing massive stars in this rather 
simplistic analysis and it is probably fair to say that the form of the CO SED observed in 
our models can be taken to be roughly representative of such active systems. 

3.5. Trends in ^^CO Line Emission at High-z 

We now consider the relative velocity integrated CO antenna temperatures for the high 
redshift models considered in B09. The physical parameters associated with these models 
are listed in Table [TJ All models represent active galaxies with high cosmic ray fiuxes and 
FUV radiation fields. Densities are also moderately high while both sub-solar and solar 
metallicities have been considered as the metallicity of high redshift systems still remains 
debated. B09 and references therein provide full justification for the choice of input param- 
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eters for these models. Once again we emphasise that these parameters are not meant to 
be representative of particular high redshift galaxies that have been observed in ^^CO but 
rather crudely represent our best guess as to the parameters that are likely to be associated 
with high redshift systems. In B09 it was found that while there is no evidence for Models 
I or II having a high-mass biased IMF, Model III may have a high-mass biased IMF if the 
magnetic pressure is sufficient to halt the collapse and formation of low-mass stars. In other 
words, the physical conditions associated with Model III are likely to be necessary but not 
sufficient for a top-heavy IMF to exist. We consider differences in CO antenna temperatures 
in these three models as a possible means of distinguishing between the two IMF scenarios 
and the differing physical conditions in these galaxies. 

Figure E] plots the relative velocity integrated CO antenna temperature for the different 
CO rotational lines for the three high-z models. Solid lines correspond to ^ 3 while 
dashed lines correspond to ^ 8. We can immediately see that the CO SEDs have very 
different forms for the three models. 

Model I has very low velocity integrated antenna temperatures compared to the other 
two models. This is because the extreme cosmic ray fluxes and FUV fields invoked in this 
model destroy most of the CO, resulting in very low line intensities. The highest J transitions 
are found to be the brightest for this model but the temperatures associated with the lines 
are always very low probably indicating that such extreme conditions may be unphysical 
even in high redshift systems. Model II and Model III both have the same densities, cosmic 
ray fluxes and FUV flelds, while Model II is at sub-solar metallicity. Model III has solar 
metallicity. The solar metallicity model is seen to produce brighter J {TA)'^^^dv than the 
sub-solar metallicity model. However, in both cases, it is the intermediate transitions that 
dominate the CO emission although the drop in the SED at high J is much sharper for the 
solar-metallicity model. At A^ ~8, the CO SED for Model II peaks at J=5 before remaining 
roughly constant whereas the CO SED for Model III drops signiflcantly at high J at the 
same A^. These trends may be used to distinguish between solar and sub-solar metallicity 
active galaxies at high redshift. Figure [5] clearly illustrates that the ratio of high-order CO 
transitions to the peak CO transition will be close to unity for physical conditions associated 
with an unbiased IMF at high redshift. However this ratio is likely to be smaller (~ 0.5) for 
physical conditions associated with a top-heavy IMF at high redshift. Unless such low CO 
J-ratios are seen in active galaxies at high redshift, an unbiased IMF cannot be ruled out. 
We caution the reader, however, that the physical conditions assumed in the models in this 
paper seem necessary but may not be sufficient to produce the IMFs associated with them 
in Table [TJ Also, in § 13.21 we inferred that sub-solar metallicity counterparts of Milky Way 
type galaxies may also produce top-heavy IMFs but in this case the ratio of the high-J CO 
lines to the peak CO line is likely to be close to unity. Larger ratios of the CO high-J to 



- 9 - 



peak transitions therefore seem to imply that the galaxy is at sub-solar metallicity and while 
such galaxies with normal levels of cosmic ray ionisation and FUV flux show evidence for 
producing top-heavy IMFs, the same galaxies at high-redshifts with elevated cosmic ray and 
FUV fluxes, may produce unbiased IMFs. The high-J to peak CO ratio therefore provides 
an useful indication of the metallicity of a galaxy and assuming it is known whether the 
galaxy is active, e.g., due to the presence of an AGN or starburst nucleus, the form of the 
IMF may be crudely inferred. 

In Figure [6] we show variations in the theoretical velocity integrated antenna temper- 
atures for the three models with increasing depth into the cloud. In all three models, the 
CO rotational lines get brighter with increasing depth within the cloud. Note that due to 
the differing number densities and metallicities of each model, the same corresponds to 
different physical distances for these models. Figure [7] illustrates the relationship between 
and size of the emitting region for our high redshift models. 

We note that the difference in the theoretical velocity integrated antenna temperatures 
of Model II (sub-solar metallicity and unbiased IMF) and Model III (solar metallicity, top- 
heavy IMF) is greater for the higher order CO transitions deep within the cloud. At A^ > 4, 
the temperatures associated with the CO(l-O) line for Model II and Model III are almost 
the same. However, the emitting region in Model II is several orders of magnitude larger 
than in Model III. The CO (7-6) or CO (9-8) transition on the other hand produces markedly 
different antenna temperatures at all depths for the three high-redshift models. 



3.6. Comparison to Observations 

As previously mentioned, ^^CO line ratios are a good way of comparing our theoretical 
predictions directly to observations. In this section, we consider line ratios for two types 
of observed extragalactic sources using where available line intensities or integrated tem- 
peratures measured through the same beam size for calculating these ratios. Once again 
we emphasise that none of our models mimic the exact physical conditions in the galaxies 
described in this section. Rather, they may crudely represent the average conditions within 
a PDR in a galaxy of that type. 



3.6.1. Nearby Starburst Galaxies - e.g. NGC253 



The nearby starburst galaxy NGC253 has been extensively mapped in -"^^ CO emission by 
various authors JSavet et allbood : [Bradford et alllioosl : Iciisten et al.lbooeh . | Bradford et al.l 
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( 120031 ) compute line intensities for various CO transitions corrected to a 15" beam. We use 
these values from their Table 3 (Column 6) to derive line intensity ratios for a few of the 
CO lines and find observed line intensity ratios of (jQ^^i^ff^ ^ 8.0 and ~ 3.3. In B09, 

NGC253 was cited as an example of a galaxy with a high cosmic-ray ionisation rate due to 
the presence of a starburst nucleus. The cosmic ray ionisation rate for this galaxy is thought 
to be ~ 800 times greater than in the Galaxy (IBradford et al.ll2003l ). The theoretical line 
intensity ratios derived from our model with a cosmic-ray ionisation rate of ( = 10~^^ s~^ 
described in § 13.31 are co\ i-ff} ^ ^'"^ co(4-3) ~ 4.1 at ~ 8. Taking the ratio of 
intensities from Table 3 of Bayet et al. (2004) we find ^qI^Io] ~ 9.3 for a beam size of 23" 

Harrison et al. 1999) in both cases and ^^l!"!! ~ 4.5 for a beam 



(iMauersberger et al. 



1996 



C70(4-3) 



size of 21.9" (iGuesten et al.l Il993l ) in both cases. The ratio of fluxes given in Table 1 of 
Glisten et al.l ( 2006 ) gives ^^l^'g^ ~ 2. Given that we have not attempted to accurately 



model the physical conditions of this particular galaxy, the agreement between the observed 
and theoretical results is remarkably good. Note that these ratios are the ratios of the 
line intensities, /, rather than th e velocity integrated t emperature plotted in the figures as 
only the intensities in Table 3 of [Bradford et al.l (120031 ) are all corrected to the same beam 
si ze. We no t e tha t the ratios of our velocity integrated temperatures agree with those in 



Bayet et al.l (120041 ) to the same level as the inte nsity ratios. Ou r pred icted line intensities 



and ratios also agree with the theoretical work of iMeijerink et al.l (120061 ) (Tables 2, 3 and 4) 
who have looked at the effect of high cosmic-ray ionisation rates and FUV radiation fields 
on the ^^CO lines 



3.6.2. High Redshift Galaxies - e.g. Cloverleaf 
The gravitationally lensed QSO, Cloverleaf (H1413 -I-117) at a redshift of ~ 2.5, was 



cited as an example of a high-redshift source in B09. iBarvainis et al.l (119971 ) report CO 
observations for this source and find the CO (4-3) line to be the strongest line in terms 
of the brightness temperature. These authors compute brightness temperature ratios of 

ggglg] = 0.83 ± 0.16, cot-f) = ^-^^ ^ co\[-3) = 0-68 ± 0-13 relative to the 

brightest line. There is therefore, a definite drop in the ratios relative to the brightest line 
when going from the low- J to the high- J transitions in this source. Figure [5] shows that for 
Model II at both ~ 3 and ~ 8 as well as Model III at ~ 3, the CO (5-4) line is 
the strongest line. For Model III at ~ 8, the CO (3-2) line is marginally brighter than 
the CO (4-3) line and is now the strongest line. We compute ratios relative to the brightest 
line for both our high redshift models II and III. For Model II which shows little evidence 
for having a high-mass biased IMF, the ratios go from ~0.7 for the low-J transitions to ~1 
for the CO (7-6) transition relative to the brighest line. For Model III which has some of the 
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physical conditions necessary for producing a top-heavy IMF, the ratios go from ~0.9 for the 
low-J transitions to ~0.7 for the CO(7-6) hne relative to the brightest line. Thus, the trends 
predicted by Model III seem to match the observations for this source better than those 
predicted by Model II. Once again, we emphasise that we have not tried, in our models, to 
match the physical conditions for this source in any way whatsoever or to try and reproduce 
the observed SED of this galaxy. However, our theoretical work does present some evidence, 
albeit speculative, that the Cloverleaf source at redshift ~2.5 has roughly solar metallicity 
and potentially a top-heavy IMF. It is expected that future observational studies will be able 
to test the validity of this prediction. 



4. Molecular Tracers of IMF at High-z 

Although the focus of the paper so far has been on the ^^CO SED associated with 
different physical conditions, the UCL_PDR code includes a network of 83 chemical species 
whose abundances are also output by the code. In this section we consider how the fractional 
abundances of some of these different molecular species that could potentially be observed 
with future surveys, change for different IMFs at high redshift. The fractional abundance 
is defined to be the abundance of the molecule relative to the total hydrogen abundance, 

= n{H) + n{H2). Once again we consider Models I, II and III of B09. In Table [2] we 
list the fractional abundances of some key species for the three models both at ^ 3 and 

~ 8 typical of the translucent gas component and dense PDR gas component of galaxies. 
In Figure [8] we also plot the variation of the molecular abundances of some key species with 
depth inside the cloud for the different high-redshift models. 

From Table [2] and Figure [HI we observe the following trends. Assuming a limit of 
detectability of 10~^^ in fractio nal abundance which is arbitrary but roughly satisfied in our 
own galaxy JSavet et al.ll2009h . all the species apart from CO are below this limit for Model 



I where the extremely high cosmic ray fluxes and FUV fields serve to destroy molecules 
effectively even deep within the cloud. Note that our qualitative results do not change if we 
assume a slightly higher or lower detectability limit. 

For Model III, the fractional abundances of all molecules considered here are above the 
detectability limit even at low Ajj where we normally do not expect complex molecules to 
be present. However, the high density in this model means that the systems considered here 
are always likely to be compact rather than diffuse and the size of the cloud is smaller than 
for any of the other high redshift models at a given Aj, as is illustrated by Figure [71 The 
high fractional abundance of CO in this model also traces a high abundance of H2 even at 
small A^. 
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While HCN and HNC are detectable in Model III at all A^, these molecules only start 
becoming detectable in Model II at ~ 2. However, the size of the emitting region in 
Model II is much larger than in Model III. At high Ay, the abundances of both HCN and 
HNC are almost the same in Models II and HI corresponding to unbiased and high-mass 
biased IMFs respectively. 

The abundance of CO and CS is found to scale linearly with the metallicity - these 
abundances are about two orders of magnitude greater for Model HI compared to Model II 
at all depths and the metallicity in Model HI is also about two orders of magnitude greater 
than in Model II. A detection of CS at a much higher level than otherwise expected in a high 
redshift active galaxy, may therefore signal that it is at roughly solar metallicity and has some 
of the physical conditions necessary for producing a top-heavy IMF. However, we note that 
the predicted CS fractional abundance may be sensitive to the assumed initial abundance 
of sulphur in our models which is taken to be 2 x 10~^ here. This fractional abundance 



reproduces roughly the obs erved abundance o f CS seen in TMC-1 (lOhishi fc Kaiful Il998l ) 



but is still highly uncertain (IRuffle et al.lll999l). These res ults for the abundance of CO and 



CS are consistent with the findings of iBayet et al.l (120091 ). These authors have studied the 



molecular tracers of various PDR dominated galaxies with different physical conditions. For 
a density of 10^ cm~^, a cosmic ray ionisation rate of 5 x 10~^^ s~^ and a FUV radiation 
field strength of x = 1.7 x 10'^, they find that CO and CS are linear tracers of the metallicity 
while CN and HCN are relatively insensitive to the metallicity. For A^ > 5, the CN and 
HCN abundance are indeed found to be roughly the same for Models II and HI. However, 
at lower A^, both CN and HCN are much more abundant in Model HI. 

The trend of HCO"^ with Ay provides another possibly interesting diagnostic for different 
IMFs at high redshift. This molecule is more abundant at low Ay {Ay < 3.5) in the sub-solar 
metallicity model which is thought to produce an unbiased IMF. However at Ay > 3.5 it 
becomes more abundant in the solar metallicity model which is thought to produce a top- 
heavy IMF. For Model HI, there is an obvious change in the chemistry at Ay ~ 3.5 and 
a sharp change in the abundance of all the molecules is seen at this Ay. It is well known 
that the attenuation of the FUV flux with distance in a photon-dominated region, produces 
a characteristic depth dependent chemical structure with well-defined chemical zones. At 
Ay ~ 3.5 in this model, we witness the transition from a cloud layer of singly ionized carbon 
produced by photo-ionization of carbon atoms by the FUV photons, to a layer where most 
of the atomic carbon becomes locked up in the stable molecule CO due to the inability of the 
FUV field to penetrate the cloud further. It should be noted that in our other high redshift 
models, which are both at sub-solar metallicities, this transition between chemical zones is 
not seen to occur within the Ay range considered in Figure [81 
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Although the variation of molecular fractional abundances with demonstrates some 
interesting differences between the various IMF models, in reality these trends may be hard 
to observe. Most observations of extragalactic sources lead to measurements of the total 
column density of different molecules integrated over the depth of the cloud. Therefore, the 
observed column density of molecules such as HCO"'" for example in two extragalactic sources 
with different IMFs may actually be very similar and therefore not a good diagnostic for the 
form of the IMF. 



5. Discussion & Conclusions 

If top-heavy IMFs arise from the suppression of low-mass star formation then the galactic 
physical conditions that inhibit star formation may be crudely defined. Physical conditions in 
which low-mass star formation is permitted are here assumed to give rise to unbiased IMFs. 
Since the physical conditions determine the chemistry, potential molecular diagnostics of 
various types of IMF can be deduced. 

The sensitivity of the ^^CO rotational spectrum to changes in gas number density, 
metallicity, and cosmic ray and far-UV fluxes are presented and used as a basis for predicting 
the relative velocity-integrated antenna temperature of the CO rotational spectrum in three 
models of high redshift active galaxies defined in B09. While active galaxies with near-solar 
metallicity (and top-heavy IMFs) produce higher antenna temperatures than those with sub- 
solar metallicity (and unbiased IMFs) the actual rotational distribution is fairly similar. In 
sub-solar metallicity counterparts of Milky Way type galaxies, the high-J to peak CO line 
ratio approaches unity. These physical conditions are thought to produce top-heavy IMFs. 
However, sub-solar metallicity active galaxies at high redshift also have high-J to peak CO 
line ratios that approach unity and these physical conditions are thought to produce unbiased 
IMFs. Solar metallicity active galaxies at high redshift have high-J to peak CO ratios that 
drop well below unity and these conditions may produce top-heavy IMFs. The form of the 
CO SED therefore clearly depends strongly on the metallicity of the galaxy and provided it is 
known whether the galaxy is active or quiescent, the form of the IMF that may be produced 
under such physical conditions, can be crudely inferred. High order CO transitions are also 
found to provide a good diagnostic for high far-UV intensity counterparts of Milky Way type 
systems which show some evidence for producing high-mass biased stellar IMFs. 

It is shown that the theoretical predictions for the ^^CO rotational spectrum can be 
matched to observed line ratios of different kinds of extragalactic sources both nearby and in 
the distant Universe. Potentially exciting new evidence is presented for the Cloverleaf QSO 
at z~ 2.5 having the physical conditions necessary for producing a top-heavy IMF. This is 
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consistent with the studies of IPavel (120081 ): IWilkins et all (120081 ) and Ivan Dokkuinl (120081 ) 



who require the IMF to be top-heavy at z > 1. 

Molecular tracers of dense gas are predicted to differ markedly between galaxies with 
top-heavy and unbiased IMFs. In particular, CO and CS differ significantly over the Ay range 
1 to 9 visual magnitudes with these molecules acting as linear tracers of the metallicity. As 
solar metallicity active galaxies at high redshift are thought to have physical conditions 
necessary for a top-heavy IMF, detection of elevated levels of CO and CS in high-redshift 
active galaxies may signal that such physical conditions exist in these sources. HCN, HNC 
and CN differ between solar and sub-solar metallicity active galaxies at high redshift for small 
visual extinctions but at large depths, these molecules are found to be independent of the 
IMF shape. Finally, HCO"'' is more abundant in active galaxies with unbiased IMFs at low 
but, at high Ay, the abundance of this molecule is greater in active galaxies with top-heavy 
IMFs. Unless the Ay of the host galaxy can be directly measured, this molecule is probably 
not a good diagnostic for the shape of the IMF as the integrated HCO"*" column density that 
would arise from our two high-redshift models with different IMFs, would probably be very 
similar. 

We conclude that chemistry may in some cases provide a simple and independent means 
of determining the approximate nature of the IMF of a high redshift galaxy. 
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Table 1: High Redshift Models: inputs to the model as well as the form of the IMF associated with these physical 
conditions and inferred in B09. Note that these physical conditions seem necessary but may not be sufficient for 
producing the IMFs associated with them. 



Model 


Density (cm ^) 


Metallicity (^q) 


CR lonisation Rate (s ^) 


FUV Flux (Draine) 


IMF 


I 


10^ 


0.05 


10-13 


1.7x10^ 


Unbiased 


II 


10^ 


0.05 


lQ-14 


1.7x10^ 


Unbiased 


III 


10^ 


1.00 


10-14 


1.7x10^ 


Top-heavy 



Table 2: Fractional Molecular Abundances for different IMF Models 



Molecule 


Model I 


Model II 


Model III 




A. ~3 
























HCN 


5.598 X 10- 


-15 


2.813 X 10- 


-14 


6.931 X 10- 


-12 


2.869 X 10' 


-9 


4.753 X 10" 


-10 


2.253 X 10' 


-9 


HNC 


6.284 X 10" 


-17 


3.498 X 10- 


-17 


1.475 X 10- 


-11 


1.770 X 10' 


-9 


4.753 X 10- 


10 


2.655 X 10" 


-9 


cs 


2.905 X 10- 


-18 


7.077 X 10- 


-18 


2.316 X 10- 


-13 


6.405 X 10- 


12 


9.496 X 10- 


10 


3.719 X 10- 


10 


CN 


3.008 X 10- 


-13 


1.153 X 10- 


-12 


2.094 X 10" 


-10 


4.474 X 10' 


-9 


4.029 X 10' 


-8 


1.758 X 10' 


-9 


HCO+ 


3.279 X 10- 


-13 


2.685 X 10- 


-13 


1.553 X 10' 


-9 


1.438 X 10' 


-8 


1.413 X 10- 


10 


7.832 X 10' 


-8 


CO 


8.037 X 10- 


-9 


1.050 X 10 


-8 


5.063 X 10' 


-6 


5.591 X 10' 


-6 


3.183 X 10" 


-5 


1.397 X 10" 


-4 



10 

2 4 6 8 

J 

upper 



Fig. 1. — Variation of relative velocity integrated CO antenna temperature per unit parsec 
with number density of hydrogen at Ajj ~3 (solid lines) and ^ 8 (dashed lines). The 
metallicity, cosmic ray ionisation rate and incident FUV radiation intensity are fixed at 
e = ^0, C = and x = 1.7. 
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123456789 



Fig. 2. — Variation of relative velocity integrated CO antenna temperatures per unit parsec 
with metallicity at ~3 (solid lines) and A^, ~ 8 (dashed lines). The number density 
of hydrogen, cosmic ray ionisation rate and incident FUV radiation intensity are fixed at 
n = lO^cm"^, ( = 10~-^'^s~^ and x = 
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7000 




upper 



Fig. 3. — Variation of relative velocity integrated CO antenna temperature per unit parsec 
with the cosmic ray ionisation rate at A^, ~3 (solid lines) and ~ 8 (dashed lines). The 
number density of hydrogen, metallicity and incident FUV radiation intensity are fixed at 
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16000 




2 4 6 8 

J 



upper 

Fig. 4. — Variation of relative velocity integrated CO antenna temperature per unit parsec 
with the FUV radiation field strength at ~3 (solid lines) and A^, ~ 8 (dashed lines). 
The number density of hydrogen, metallicity and cosmic ray ionisation rate are fixed at 
n = lO^cm"^, ^ = ^0 and ( = 1 x 10~^'^s"^. 
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Model I 




Fig. 5. — Variation of theoretical velocity integrated CO antenna temperatures for high 
redshift models (Table [T]) at A^, ~3 (solid lines) and A^, ~ 8 (dashed lines). 
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Fig. 6. — Variation of theoretical velocity integrated CO antenna temperatures with for 
high redshift models (Tabled]). 
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Fig. 7. — Relationship between the size of the emitting region and the visual extinction, Ay 
for our three high redshift models. The dotted lines correspond to the two AyS at which we 
consider our outputs. 




Fig. 8. — Variation of fractional abundances of different molecular species with depth, 
for high-redshift models (Table [1]). The solid black line indicates the assumed limit of 
detectability of lO^^^ 
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